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ABSTRACT 

In the Halo Model, galaxies are hosted by dark matter halos, while the halos them¬ 
selves are biased tracers of the underlying matter distribution. Measurements of galaxy 
correlation functions include contributions both from galaxies in different halos, and 
from galaxies in the same halo (the so-called 1-halo terms). We show that, for highly 
biased tracers, the 1-halo term of the power spectrum obeys a steep scaling relation in 
terms of bias. We also show that the 1-halo term of the trispectrum has a steep scaling 
with bias. The steepness of these scaling relations is such that the 1-halo terms can 
become key contributions to the n-point correlation functions, even at large scales. We 
interpret these results through analytical arguments and semi-analytical calculations 
in terms of the statistical properties of halos. 
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1 INTRODUCTION 


Galaxy surveys l|York et al.||2000| |Cole et al. 

2005| Abbott 

et al.|2005l|Scoville et al. 2007l|PAN-STARRs| 

Dawson et al. 

2012 

Blake et al.|2011 

Anderson et al.|2012 

2014 

1 are not 


just tools for constraining cosmological parameters: they re¬ 
veal the 3-dimensional spatial web of visible structures, the 
time evolution of these structures, and the history of galaxy 
evolution. Next-generation surveys are aimed at answering 
a variety of open astrophysical and cosmological questions, 
by collecting vast amounts of data at low, intermediate, and 


high redshifts (Schlegel et al.|2009t Benitez et al.|2009| 

2014 

Ellis et al. 

2012||Dawson et al.|2015| |Levl et al.|2013| 

Ib^ 

et al.|2009 

1 . 


For the study of how gravity and the Universe’s back¬ 
ground expansion affect the growth of structure, the stan¬ 
dard statistical tools are the n-point correlation functions 
and their Fourier transforms, the polyspectra. If the dis¬ 
tribution of matter were perfectly Gaussian, the two-point 
correlation function (2pCF) or, equivalently, the matter 
power spectrum, would contain all the statistical informa¬ 
tion. However, the primordial fluctuations are believed to 
be very nearly, but not perfectly, Gaussian. Information 
about the processes that generated these primordial per¬ 
turbations is encoded in higher order moments, e.g. the bis¬ 
pectrum ( Maldacena|2003' \. Moreover, when subject to non¬ 
linear time evolution, even a perfectly Gaussian Initial held 
develops nontrivial higher-order moments (skewness, kurto- 
sis, etc.) whereto statistical information propagates. At late 
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stage of non-linear gravitational evolution, information even 
leaks out of the hierarchy of moments, as the density Held 
becomes approximately lognormal (Garron|2011 2012 Car 


ron fc Neyrlnck|2012| Carron|[2014 1. 

Adding to these complications is the fact that we do not 
directly observe the total matter distribution, but only its 
visible component — which accounts to ~ 20% of the total 
matter (see, e.g. Planck Collaboration|20l5 \ and is affected 
by non-gravitational effects such as the physics and feedback 
of baryons, radiation pressure, etc. Hence, the 2pCF of the 
distribution of visible matter cannot possibly tell the full 
story and we must treat galaxies, quasars, Ly-a systems, 
etc., as unfaithful (and biased) tracers of the underlying dark 
matter (DM) distribution. 

The relation between tracers of large-scale structure 
(LSS) and the DM distribution is partially provided by the 
Halo Model ( Cooray fc Sheth||2002 1. The DM halos — and 
not the DM particles — then become the fundamental ob¬ 
jects. In particular, the correlation functions of the DM halos 
are related to the correlation functions of the DM particles 
by the halo abundance, bias and proHle, such that more 
massive halos are less abundant, are more highly biased and 
have less concentrated proHles. 

The statistics of how galaxies populate DM halos is pro¬ 
vided by the so-called Halo Occupation Distribution (HOD) 


( Martinez fc Saar|2001| Berlind fc Weinberg|2002 Kravtsov 


et al. 2004 Zheng et al. 20051, which specify how many 


observable galaxies inhabit a halo as a function of its mass. 
HOD parameters can be calibrated by measurements of spe- 
cihc abundance ratios, and by Htting the observed correla¬ 
tion functions. 
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The relationship between the statistics of the halo den¬ 
sity field and that of galaxies is often non-trivial. For ex¬ 
ample, the n-polnt statistics of galaxies get contributions 
from the A^-halo term (when each galaxy occupies a differ¬ 
ent halo), from the (A'’ — l)-halo term (when two galaxies 
occupy the same halo, and the others lie in different ha¬ 
los), etc., all the way down to the 1-halo term (when all the 
N galaxies occupy the same halo). Furthermore, there are 
contributions to the A^-point statistics from A'^ — 1 types of 
shot-noise terms: e.g., the {N — l)-halo - (A'^ — l)-galaxy 
term (when the correlation function hits twice a galaxy in 
a given halo, and then only once galaxies in different halos) 
down to the 1-halo - 1-galaxy term (when the correlation 
function hits N times the same galaxy). Hence, when we 
measure the A^-point function of galaxies, we are in fact 
measuring an ad-mixture of all the A^'-point functions of 
halos {N' = 1,2,..., AT). 

The galaxy 2pCF includes both a 2-halo term (related 
to galaxies in two distinct halos) and a 1-halo term (account¬ 
ing for galaxies within the same halo). Clearly, the 2-halo 
term Is most important on large scales, within the linear 
regime, and bears the imprint of the large-scale matter dis¬ 
tribution, whereas the 1-halo term dominates on small, non¬ 
linear scales, and reflects the matter distribution inside halos 
(the halo density profile). Even though these two regimes are 
connected by time-evolution, we may describe the galaxy 
power spectrum as the superposition of two independent 
scaling laws; that of the linear DM power spectrum, which 
dominates on large scales, and that coming from the halo 
profiles, which dominates on small scales {k > Ih Mpc“^). 

In the large-scale limit the 1-halo term contributes a 
constant to the galaxy power spectrum. This constant may 
be poorly known, since galaxy surveys designed for cosmo¬ 
logical studies are often insufficiently complete (or have poor 
redshift accuracy, in case of photometric redshifts) to deter¬ 
mine precisely the HOD parameters for a type of galaxy, 
and for the survey’s mean redshift. This constant represents 
a noise that must be subtracted from the power spectrum 
— just as it happens with shot noise, which has a com¬ 
pletely different origin but is also a constant that must be 
subtracted. This feature of the 1-halo term on large scales 
also comes into play in higher-order statistics: the bispec¬ 
trum gets a constant contribution from the 1-halo term of 
the 3-polnt function; the trispectrum gets a constant contri¬ 
bution from the 1-halo term of the 4-point function; and so 
on. 

In this paper we study the behavior of the 1-halo terms 
as a function of the galaxy bias, where these quantities are 
connected via their mutual dependence on HOD parameters. 
The same goes for the other observable quantities, such as 
the mean number density of galaxies, the 2-halo term of the 
power spectrum, etc. Hence, the HOD provides an internal 
(but unobservable) parameter space that we can use to vary 
the observable quantities. In particular, we use the galaxy 
bias to parametrize the 1-halo terms both because it is more 
readily available in observations, and also because we are 
interested in identifying universal behaviors, regardless of 
the details of the HODs. 

We show that, for highly biased tracers (6g > 3), the 
1-halo term of the 2pCF obeys a scaling relation in terms of 
bias, growing as ~ which is much faster than the 

scaling of the 2-halo term (P^'* = b^Pm, where Pm is the 


matter power spectrum). For highly biased galaxies, the ef¬ 
fective shot noise contributed by the 1-halo term can become 
at least comparable to the Poisson shot noise, significantly 
lowering the signal-to-noise ratio for measurements of the 
power spectrum, baryon acoustic oscillations, etc. In some 
cases the 1-halo term can even surpass by a large factor the 
shot noise, as is the case, e.g., for the angular power spec¬ 
trum of the cosmic infrared background on the angular scales 
probed by Herschel ( Thacker et al.|[2013 l. Furthermore, we 
show that the 1-halo term for the trispectrum also grows 
very fast — typically, like (pi'*)^ — and should, therefore, 
contribute an important source of noise for the power spec¬ 
trum covariance in the limit of high bias. 

When employing a particular cosmological model be¬ 
low, we used a standard flat ACDM scenario, with Qm ~ 
0.26, ris = 0.96 and ag = 0.78. 


2 FORMALISM AND ANALYTICAL 
APPROXIMATIONS 

2.1 The Halo Model 

Over time, gravity enhances the density contrast field by 
attracting matter towards the density peaks, and by cre¬ 
ating voids where the density was initially below-average. 
The Halo Model ( Cooray fc Sheth|2002 I describes how this 
process depends on the mass of the collapsed structures by 
determining, e.g., the abundances of the DM halos — i.e., 
the mass function, dfih/d logM. 

Halos are ultimately determined by the peaks of the 
initial density field, and according to the theory of peak 
statistics ( Bardeen et al.[1986 Mo fc White|199'6 1, the main 
driver of the abundance of peaks as a function of mass is 
the mass variance within a certain comoving radius R. The 
variance of the linear density field inside a spherical top-hat 
region of radius R is determined by: 


a^{R) = ^ f dkk^Pm{k)W{kR) 


( 1 ) 


= / 


d\nk/Fi,{k)W{kR) 


where Pm{k) is the linear matter power spectrum and 
W{x) = \Zj\{x)/x^ is the window function for a spheri¬ 
cal top-hat region. The mass contained within radius R at 
the mean background matter density today, pm, is M{R) = 
4nR^ pm/3. The peak height is defined as v{M) = Sc/cr{M), 
where Sc is the linearly extrapolated critical density contrast 
for spherical collapse. 


Due to the statistics of density peaks (Bardeen et al. 


19861, all mass functions exhibit an exponential dependence 


on the peak height v{M). In fact, we have: 


duh p m diner 


d\nM M d\nM 


/H 


( 2 ) 


where, up to some model-dependent factors and coefficients, 
f{v) ~ 

Another crucial ingredient of the Halo Model is the 
halo bias. Since we would like to substitute the true mat¬ 
ter density contrast 5m by the contrast of halo counts, 
Sh = nh/rih — a relationship between the two must be 
established. We can write this in terms of a local ansatz 
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such as (Fry fc Gaztana^|1993 1: 


Sh = hh5m + - cjf) + ■ 


(3) 


In this paper we will only consider the first term in this 
relation — the higher-order terms can also become impor¬ 
tant precisely in the limit that we are investigating (tracers 
with high bias), but we leave this key issue for future in¬ 
vestigations. Typically, the halo bias is a smooth power-law 
function of peak height, such that more massive halos cor¬ 
respond to higher (and rarer) peaks and have higher values 
of the halo bias. 

Our analytical calculations were performed using both 
the Press-Schechter ( Press fc Schecht^|1974|) (PS) as well 
as the Sheth-Tormen fSheth k. Tormen||1999[ [Sheth, Mo| 
fc Tormeii] 19991 (ST) prescriptions for the mass function 
and halo bias. The semi-analytical calculations of Sec. |3.3 


were performed using the Tinker mass function ([Tinker et al. 


20081 and halo bias (Tinker et al. 20101. The main results 


are very similar on all cases. 


2.2 Halo Profile 

Related to the Halo Model, but still a slightly orthogonal re¬ 
sult which depends more strongly on the non-linear regime 
of structure formation, is the shape of the halo density pro¬ 
file. Although our results are completely insensitive to the 
fine details of these profiles, in our semi-analytical calcula¬ 
tions we employ the standard results of [Navarro, Frenk fcj 
White (19971. In that case, the density profile for a halo of 


mass M is given by: 

p(r|M) ~ u(r\M) 




(4) 


where = rs{M) is the characteristic scale (“knee”) for a 
halo of mass M. The mass-averaged halo profile in Fourier 
space is computed as: 

drr —,-p(r|M) . (5) 

kr 


<t\M) = I / 


The important feature for our purpose is that u{k\M) —> 1 
when k < Ih Mpc“^ for the range of masses that we are 
interested in. This property follows simply from the fact 
that u{k —>■ 0) = (1/M) f d^xp = 1, and in the range k < 
1 h Mpc“^ this approximation holds for the mass scales we 
are interested in. Hence our calculations are insensitive to 
the precise adopted profile, and importantly to whether or 
not the galaxy count profile is identical to the DM density 
profile. Note indeed that a discrepancy between these two 
profiles may have been observed in the innermost regions of 
the DM halos, where the galaxy profiles may be significantly 
steeper than the NFW profil e ([Watson et al.|20Toj [Kayo fc| 
Qguri|20l2 Piscionere et al. [[20141. 


2.3 Halo Occupation Distribution 

It has been known for a long time that there must be a direct 
relation between the distribution of galaxies and that of the 


underlying dark matter halos ( 

White fc Frenk[[1991 

Kauff- 

mann. White fc Guideroni||1993 Navarro, Frenk fc White 

19951 

Mo fc White 1996| |Kautfmann et al.||1999| [Springel 

et al. 

2005 

I. Halo Occupation Distribution (HOD) models 


provide the probability distribution function P{N\M) for a 
certain number (N) of galaxies to occupy halos of a given 
mass (M) (Ma & Fry 2000 [Seljak 2000 Cooray & Sheth 


2002| [Martinez fc Saar[[200i| [Her ind 6i Weinberg[ [200^ 
Zheng et al. 2005 I. Although halo mass is not the only factor 


which determines the number of galaxies (Zentner, Hearin & 
van den BoschpOld I, HODs have been extremely useful to 


interpret measurements of the clustering of different types 


of tracers, from galaxies (Zheng, Coil fc Zehavi|2007 Zheng 

et al. 

20091 to quasars 

(Porciani, Magliocchetti, fc Norberg 

2004 

Shen et al.||2007| 

Wake et al. 2008| Shen et al.||2010| 

Kayo fc Oguri 2012 1 

Lichardson et al. 2012 1 — however. 


in the latter case the simplest HODs may be inadequate to 
capture the complex interactions between quasars and their 
environments ( Shen et al.|2013 Chatterjee et al. ]M^[^ 


fc Safarzadeh[[20T5 1, and additional parameters such as as¬ 

sembly bias should be included. Often, instead of P{N\M), 
what is provided are the momenta of the HOD, such as 
N{M) = {N)m, {N{N — 1))m, etc. The brackets define av¬ 
erages over halos of the same mass, and the HOD can be 
defined in terms of these momenta. For brevity, we will drop 
the subscript M from now on. 

It is clear that, for very massive halos, the number of 
galaxies should scale proportionally to the halo mass, but 
as we approach the low-mass end the situation can be more 
nuanced. According to the hierarchical scenario of structure 
formation, a galaxy can either form inside its original halo, 
or join after formation an already existing and populated 
halo. Such a dichotomy is also seen in numerical simulations 
( [Kravtsov et al.[[2004[ ) and leads to the distinction between 
“central” galaxies, of which there can be only one per halo, 
and possibly numerous “satellite” galaxies. A popular func¬ 
tional form for the number of central and satellite galaxies 


is (see, e.g., Zheng et al. (2005l): 


Vx 1 / Me - M 

(JV.) . A'.-5Exfc(^j- 

{Ns) = Ns=NsXNs , 

where: 


,x ex., M — KqMc 

Ns = 9{M - KgMs) ' ® 


Ml 


( 6 ) 

(7) 

( 8 ) 


As denoted by Eq. ([^, the existence of satellites is con¬ 
ditional on the existence of a central galaxy. Typical val 
ues for the HOD parameters are Me — 10 


13-5 /i-i 


Ml ~ h ^ Mq, a ~ 0.9 — 1.0, Kg ~ 1.1, o-g ~ 1 (Zheng 


et al. 2005 \ — although, in the case of quasars, especially 


at high redhifts, some parameters can deviate significantly 
from these values ( [Chatterjee et al.|2013| . 

Besides the mean numbers (or richness), we must also 
specify the higher-order momenta of P{N, M). If we are only 
interested in the 2-halo and in the 1-halo terms, then all we 
need are the expectation values (N/), {N^Ns) and {Ns). The 
model separating central and satellite galaxies naturally pro¬ 
vides these momenta. By definition, the central galaxy either 
exists {Nc = 1) or does not exist {Nc = 0) inside a halo so 
{Ne{Nc — 1)) = 0 or equivalently {N^) = Nc- Regarding 
the cross-correlation between central satellite galaxies, no¬ 
tice that satellites can only exist if there is already at least 
one central galaxy, so {N^Ns) = Ns- As for the satellites, we 
can assume a simple Poisson distribution, which means, in 
particular, that {N^) = Ns{Ns + 1). 
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2.4 Combining the Halo Model and the HOD 


The Halo Model allows us to compute several quantities of 
interest from these ingredients. The mean galaxy number 
density is: 

n, = JdlnM^^xNiM), (9) 

where N = Nc + Ns. The galaxy bias is given by: 

feg(fc) = ^ J X N{M)b{M)u{k\M) , (10) 

where recall that u{k\M) —>■ 1 for k <lh Mpc“^. In terms 
of the galaxy bias, the two-halo galaxy power spectrum is 
given by: 

P^^{k)^bl{k)P^{k) . ( 11 ) 


The 1-halo term, on the other hand, is given by the correla¬ 
tion of two different galaxies in the same halo: 


P^'*(A:) 


— / dluM 

nl J 


= — / dluM 

I 


dnh 

d\nM 

duh 

d\nM 


X {N{N -1)) |u(fc|M)r 

( 12 ) 


xNAM)[2Ns{M) + N^s{M)] \u{k\M)\^ 


where the term inside square brackets in the third line is the 
intra-halo number variance, given the assumptions outlined 
above. As usual, the 1-halo term of the power spectrum does 
not include the contribution arising from self-correlations of 
galaxies with themselves (i.e., shot noise): in fact, when there 
is a single galaxy in a halo, that is, by definition, the central 
galaxy, so the number of satellites is zercQ The measured 
galaxy power spectrum is therefore: 

Pg = P^^ + -b Ps , (13) 


where Ps takes into account the shot noise power spectrum 
— which, under the assumption of Poissonian statistics for 
the galaxy counts, is given by Ps = 1/hg. 

Similar arguments can also be applied to higher-order 
correlations. The trispectrum r(A:i, A: 2 , fcs, ^ 4 ), i.e. the 4- 
point function in Fourier space, is of particular interest as 
it determines the covariance of the power spectrum, in its 
limit T{k,—k,k',—k'). The 1-halo term of this part of the 
trispectrum is given by: 


Tg^{k,-k,k',-k') 


4 f 

nj J dlnM 


n 


(14) 


x{N{N - 1){N - 2){N - 3)) 
\u{k'\M)f 
dfih 

I Cl ill ivi - 

nt 


X \u{k\M)Y 
b’j f dlnM-,, _ 

n'jj J dlnM 


xnmn! + Nt] 

X \u{k\M)Y |u(fc'|M)|'' , 

where on the second line of the equation above we used 
the same assumptions about the statistics of central and 


^ Notice that the one-halo terms vanish in the absence of satellite 
galaxies. However, due to the Poissonian nature of the HOD, it is 
possible for a halo to have satellites even if (N) = 1 for the mass 
of that halo. 



Figure 1. Power-law index as a function of M. Solid (black) 
line: n = —3(1 + dlnF/din M). Dashed (red) line: np = 
din P(k) / dink, evaluated at fc = kp = 1/R{M). The wiggles 
seen in np are caused by the BAOs. 


satellite galaxies that were used to obtain the expression in 
the second line of Eq. (12 1 . 

For our purposes, we will consider scales larger than the 
size of the largest halos, so we take u{k\M) —>■ 1 in all our 
expressions from now on. We have checked that this is an 
excellent approximation for k < Ih Mpc“^. 


2.5 Analytical model for the peak height 


In some of the following subsections we will perform analyt¬ 
ical computations of several quantities of interest in the PS 


(Press & Schechter 19741 and ST (Sheth & Tormen 1999 


Sheth, Mo & Tormen 19991 formalisms. In order to carry 


out those calculations we need an analytical approximation 
for the peak height in terms of the halo mass. 

The variance of the linear density field inside a tophat 
spherical region of radius R was given in Eq. 0 - Since 
the spherical tophat window function has the features that 
W{0) = 1 and W —>■ 0 for x S> 1, with a full width at 
half maximum of approximately Xfwhm ~ 1, it is fair to 
approximate the variance as: 


A^(fc = kp) (X k\ Pm{kp) , 


(15) 


where kp = 1/R. For self-similar models in which Pm.{k) oc 
fc", we have F{R) oc oc . In terms of the mass 

contained inside radius R at the mean background density, 
M oc R^, we have: 

(j^(M) oc . (16) 


Even though LSS in a standard ACDM Universe is not de¬ 
scribed by a self-similar model, it will be useful to consider 
this case as it will allow us to obtain interesting analytical 
expressions for fixed n. 

The peak height, normalized to 1 at M = M*, is given 
by: 


u = (M/M.)("+^>/® , (17) 

where M* is typically ~ 2. x 10^^ h~^ Mq in the ACDM 
models. In particular, within this approximation we have 
dlna~^/din M = (n-|-3)/6. 

We check the approximation of Eq. (16 1 in Fig.[^ where 
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we plot the power-law index n = —3(1 -I- d\ua^/d\nM) 
(solid, black line) together with np = d\nP{k)/d\nk, eval¬ 
uated at fc = fcp = 1/R{M) (dashed, red line). The slope of 
the power spectrum, np, shows the wiggles from the BAOs. 
The power index of the peak height, on the other hand, is an 
average over several different scales, hence it is only sensitive 
to the mean slope of the power spectrum. It is clear that the 
two are closely related, and that the approximation of Eq. 

(171 holds quite well for n between —2 and —1 in the mass 


range considered. In Sections |3.3| and |3.4| we do not use this 
approximation anymore, and instead compute a{M) from 
the power spectrum in a ACDM model. 


2.6 Mass functions and halo bias 


The simplest case is that of the Press-Schechter (PS) for¬ 


malism (Press & Schechter 19741. It provides closed-form 


expressions for the mass function and halo bias: 


fps{y) = 


bps{i^) = 1 + 


^ r/exp [—1/^/2] , 

TT 


(18) 

(19) 


These formulas are in poor agreement with the data and N- 
body simulations, however, when used in conjunction with 
an extremely simple HOD, they yield simple, straightfor¬ 
ward analytical calculations whose results convey the basic 
message of this paper — see Section |3.1| 

A better fit to simulations and observations is given by 
the ST mass function and halo bias ( |Sheth fc Tormen|1999 
Sheth, Mo & Tormen 19991: 


fsT = — [^ + {aP^) 

V TT '■ ■’ 

2p 


( 20 ) 

bsT = 1 + 3 . . , ( 21 ) 

5c S4l-e{ap^)p] 

where A ~ 0.322, a ~ 0.71, and p ~ 0.3. The ST framework 
gives a more accurate description compared to PS, while 
still allowing for fully analytical calculations. In Section [3.2| 
we use the ST formulas and a slightly more realistic HOD 
compared with the calculation in the PS case — yet the 
main results of that Section are basically unchanged. 

Finally, we also consider the expressions found by |Tin-| 
|ker et aL] ( |2008[[2010p , which were calibrated from numerical 
simulations. The mass function ( Tinker et al. [20(18 1 and bias 
(Tinker et al.|2010| are given in this case by: 


/(o-) = 0.186 X 




bh{v) = l-Aj 


-b 


-b , (22) 


-b -b , (23) 


where At = 1 -b 0.24j/ exp[—(4/i/)'*] (with y — logj^Q A, 
where we choose A = 200), at = 0.44 1 / — 0.88, Bt = 0.183, 
bT = 1.5, Ct = 0.019 -b 0.1071/ + 0.19 exp[-(4/t/)‘^], and 
ct = 2.4. We will employ this mass function in our semi- 
analytical calculations, assuming now a realistic HOD — 
see Section [ 3 .3| As we shall see shortly, the results are qual¬ 
itatively identical to those of Sections |3.1| and |3.2| which 
were found by means of analytical calculations. 


3 APPLICATIONS 

3.1 Press-Schechter mass function and a simple 
HOD 


We begin assuming an extremely simplified HOD, which 
should hold in an approximate sense for sufficiently high 


halo masses (see, e.g., Porciani, Magliocchetti, & Norberg 
(2004])): 


N{M) = 


M 

Ml 


e{M - Ml) 


(24) 


where 9{x) is the Heaviside step-function. In this simple 
HOD we take the cut-off mass to be equal to the mass scale 
Ml. This HOD also assumes that all galaxies are satellites. 
In the final subsections we recover the full description in 
terms of W and Ns, and show that the central galaxies are 
unimportant in the limit we are interested in {bg > 3). Ex¬ 
cept for the low-mass limit, the halo richness should scale 
roughly proportional to its mass, so a « 1. We will as¬ 
sume for the moment that Mi also defines the threshold 
for finding galaxies in halos - i.e., N — 0 for M < Mi. 
This approximation will be improved in the next subsec¬ 
tion, where we carry out the same calculations as here, but 
using the Sheth-Tormen formalism. As we will see, this does 
not change signihcantly our main results. 

We start by computing the number density of halos 
which host at least one galaxy in our simple HOD, Eq. ( |24[ ). 
Since the number of galaxies in each halo follows a Poisson 
distribution, this is given by: 

nu/g ^ J d In M[l - exp(-iV)] (25) 


Another definition, which will become more useful later on, 
is the number of halos that could contain galaxies: 


nh,, 


I 

I 


dlnM 


duf, 


All 
2 p. 
TT M, 


dlnM 

Pmdlna~^ [2 r 2/01 
dlnM — z/exp — z/ /2 

M dlnM V TT 


(26) 


f 


dv V 


-e/{n+3)^-N/2 


where we have used the PS mass function and the ap¬ 
proximations M = as well as the definition 

ui = (Mi/M*)("+3)/®. 


For high Mi the integral in Eq. (251 is dominated by the ex¬ 


ponential behavior of the mass function, and we can replace 
the N in the exponent by the mean number of galaxies in the 
halos just above the cut-off mass scale, N{M) —>■ N{Mi) = 
-^min- Hence, the actual number of halos containing galaxies 
can be approximated by n/^/g Ri (1 — e“^“‘") x n^.g. For the 


HOD of Eq. (241 this minimum mean number of galaxies is 
Nmin = 1, so fih/g 


0.63 nh,g. 

With the variable change x = /2 we obtain: 


(27) 

where Ao = —1/2 — 3/(n -(- 3), xi = z^i /2, and r(K, x) is the 
(upper) incomplete Gamma function of order k. Typically, 
—2 < n < — 1 for halos at the scales of interest, which means 
that -7/2 < Ao < -2. 
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The incomplete Gamma function is related to the simple 
Gamma function by r(K) = r(K, a; = 0), and has asymptotic 
limits given by: 


limr(K, a;) —>■ r(K) — a;" 

a;—^0 


K 1 + K 


+ 0{x^) 


and 


lim r(K, a;) —>■ e ^ x'^ 


1 + + 0{x-^) 

X 


(28) 


■ (29) 


Notice, in particular, that lima;_).oo r(l+ k, a:)/r(K, a:) —>■ 1 + 
a; + Cl(a:“^). It is also interesting to note that, for 0 < «: < 1, 
the asymptotic expression of Eq. ( |29[ ) is remarkably accurate 
down to a; ~ 1. 

Hence, for the ranges of interest for n, in which Aq is 
negative, the number density of the halos that host galaxies 
should diverge in the limit xi —> 0 — i.e., when Mi <C M,. 
Indeed, the number of halos of arbitrarily small masses is 
arbitrarily large, unless we specify a smoot hing scale Rf, 
in which case it asymptotes to nh,g oc (Bardeen et al. 

TMel. 


Similarly as was done above, we can compute analyt¬ 
ically the quantities defined in Section |2.4| For the mean 
number density of galaxies we obtain: 

duh 


I 


dlnM 


Ml 

2 Pm 
TV M* 


dlnM 


N{M) 


(30) 


AR 

Wi 


2-'M(l-bAi,a:i) 


where we have used the same definitions as above, with the 
difference that now the index is Ai = Ao -I- 3Q/(n -|- 3) = 
—1/2 -I- 3(a — l)/(n -I- 3). Since a ~ 1 and — 2 < n < — 1, we 
have Al ~ —1/2. Notice that for the case a = 1, Ai = —1/2 
and in the limit xi —>■ 0, we have fig = pm/Mi = pm{N)/M 


Interestingly, using Eq. (291 we find that in the high 
mass limit {xi —>■ oo), fig = fih,g - 


in that case the 


number of halos that could host a galaxy is equal to the 
mean number of galaxies. This is a consequence of the simple 


HOD, Eq. (241, which takes the cut-off mass to be identical 


to the mass scale Mi. Since the galaxy bias increases with 
Ml, for high values of this mass the number of halos above 
the cut-off is exponentially suppressed, and only the least 
massive halos are populated with galaxies. In this case each 
halo ends up hosting only one galaxy — or, more accurately, 
because of Poisson statistics, about 63% of halos contain 
only one galaxy, and the rest contain two or more galaxies. 

We can calculate the galaxy bias in the same fashion. 


using the PS halo bias of Eq. (191: 

-1 , ^ c-i r(2 -I- Al, xi) 

^ r(l + Ai,a:i) 


fag = 1 - 


(31) 


In the *1 —>■ 0 limit we can use the property of the Gamma 
function r(2 -|- Ai) = (1 -I- Ai)r(l -|- Ai), which leads to bg ~ 
1 — -I- 2dJ^(l -I- Al). We also note that taking a = 1 

leads to Ai = —1/2 and bg = 1, simply reflecting the halo 
bias consistency relation. On the other hand, in the limit of 
very large threshold masses {xi ^ 1) we obtain &g ~ 1 -|- 
4-2(5“^ xi. Hence, in order to increase bias, it is sufficient 
that xi 3> 1 — however, this is not a necessary condition: 
one could also fix the cut-off mass scale and decrease n, or 


increase a. 


A similar calculation as the one performed in Eq. (301 


leads to an expression for the 1-halo term of the galaxy 
power spectrum: 


pAn _ 


1 I ,, ., duh 

/ dlnM ——— 

J dlnM 




— 2 
TT M* 


2 A 1 -A 2 


r(l + A 2 , ail) 

[r(l-kAi,a;i)12 


(32) 


where on the first line the term {N{N — 1)) reduces to 
since we only have satellite galaxies with a Poisson distribu¬ 
tion, and on the second line we have substituted the expres¬ 
sion for hg and A 2 = Ao -I- 6a/{n + 3). Using the expression 
for the number density of halos hosting these galaxies, Eq. 
(261, we obtain: 

pih _ r(l + Aq, ail) r(l -I- A 2 , 3 : 1 ) ^ ^ 

[r(l-I- Ai,a;i)]2 Uh. 

Notice that the three indices, Ao (that appeared in fih,g), 
Al (which appeared in hg and bg), and A 2 , are related by 
Ao = 2Ai — A 2 . By substituting the asymptotic expression 
for r(K, z) in the limit 2 ; —>■ 00 one can verify that the pref- 


(33) 


actor appearing Eq. (33l is r(l-|- Ao, a:i)r(l-I- Ai, a;i)/[r(l-|- 
Ai,a;i)]^ ~ 1. Hence, in the limit of high threshold mass, 
« l/n;.,g . 

In the opposite limit, of low threshold masses, the pref¬ 
actor can become quite large — but then so does fih,g be¬ 
come large. In general, this case corresponds to tracers with 
low biases. In that limit, it is convenient to revert to the 


original expression, Eq. (321, and write instead: 


pAtl _ 


2^2-Al 

\“ r(i-bA 2 ,xi) 

llg 


I r(l-fAi,xi) 

2^2-Al 

r(i-kA 2 ) 

Ug 

\Mi, 

1 r(i-f Al) ■ 


(34) 


Let us suppose that we can change the HOD parameters 

fixed. 
-1 


while maintaining the mean number of galaxies, rig. 


In the limit of low-mass threshold the bias is 6g « 1 — <5“^ -|- 
-I- 6(a — l)/(n -I- 3)] — i.e., higher values of a and/or 
lower values of n correspond to higher biases. Since A 2 = 
Al -I- 3a/{n + 3), the ratio of Gamma functions in Eq. ( |34[ ) 
can be regarde d a s a steep function of bias. 

From Eq. ( 33 1 we see that , which is proportional to 
l/h-ii,g, plays the role of a halo shot noise term. For highly 
biased tracers (whose mass thresholds are relatively high), 
the galaxies end up in just a few halos, so that the halo 
shot noise term becomes an important part of the power 
spectrum and its covariance. This means, in particular, that 
the barrier for measuring the (2-halo) power spectrum with 
highly biased tracers may not be just the shot noise of that 
tracer, but also an additional halo shot noise coming from 
the 1-halo term. 

In the limit of high bias (and high mass thresholds) we 
can neglect the subdominant dependence in the prefactor of 
l-lhh,g in Eq. (331, and express the contribution of the halo 


shot noise term as a function of bias. Using bg ~ 1 -|- 5^, ^ -|- 


26Fxi 2 ( 5 ; 
expressed as: 


xi, the mean number density of halos can be 


nh,. 


hence: 


n + 9 

2(»» + 3) -XI 


^1+9 

,2(71 + 3) XI 


71+9 

bg + ( 


,~^g ^c/2 


ti+9 

2(ri+3) bg 5c/2 


(35) 


(36) 
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Figure 2. Scaling of the 1-halo term in the PS model. The exact 
formula, Eq. ( |32[ l, is denoted by the thick solid lines (upper line: 
n = —2; lower line: n = —1), while the dashed lines correspond 
to the approximation of Eq. | |36| l. The dotted lines are the power 
laws bg (for the case n = —1) and 6® (for the case n = —2). 


For intermediate values of the galaxy bias the 1-halo term 
is well approximated by a power-law. In Fig. we show the 
exact formula, Eq. ( |32[ |, for the cases n = — 1 (lower solid 
line) and n = —2 (upper solid line), in arbitrary units. The 


first approximation in the middle of Eq. (361 is plotted as 
the dashed lines for the two cases, where we used xi « (1 — 
Sc Sc bg)/2 . One can see that the approximation becomes 
better for higher values of the bias. Also plotted (dotted 
lines) are the power laws bg (for the case n = — 1) and 
(for the case n = —2). The simple power laws shown in Fig. 
[^are a good approximation in the interval 1.5 < bg < 3. 

Since the 2-halo term scales as = bg Pm , but the 
1 -halo term grows much faster with bias, the latter com¬ 
ponent should become increasingly important for highly- 
biased tracers. In fact, this already happens at small scales 
(fc > lb Mpc“^) even for galaxies with relatively low bi¬ 
ases. If we select tracers with increasing values of the bias 
(e.g. quasars), the 1-halo term will become more important, 
even in the large-scale limit, acting effectively as a type of 
“halo shot noise”. However, in contrast to the usual situation 
where shot noise can be beaten down by observing a larger 
number of galaxies, when the tracers are very highly biased 
this halo shot noise cannot be lowered, and a limiting factor 
for measuring the power spectrum is the finite number of 
halos, and not only the number of galaxies in the survey. 

It is easy to see that this argument also applies to 
higher-order correlation functions. The same type of inte¬ 
gral computed above appears also in the 1-halo term of the 
trispectrum, Eq. (141, and taking k —>■ 0 and k' —> 0 leads 
to: 


/^Ifl _ 


4 f dlnM-^[4iV3 + iV"i 
n* J dlnM^ 


1 / 2 Pm 

_ X 1 / _ - _ 

ni \ TT M* 


3a 




M, 


4a 


+ 1 — ) 2"^r(l + A4,a:i) 


(37) 


where Ai = Ao -I- i x 3a/(n -|- 3). In the high-mass, high-bias 


limit we obtain that: 



Thus will become a dominant part of the power spec¬ 
trum covariance matrix in the high bias limit. 

The same calculation can be employed generally for the 
1-halo term of the N -|- 1-th order polyspectrum in the 
high bias limit, showing that it grows with the scaling 

N TI.+9 -- 

1^2 n-\-3 ^Nbg 5c/2 ^ 

Hence, we conclude that for highly biased tracers not 
only the power spectrum, but also the higher-order statis¬ 
tics, are increasingly affected by intra-halo statistics, and 
may become effectively limited not only by the counts of 
the tracers, but by the counts of the halos as well. 


3.2 Sheth-Tormen mass function and a simple 

HOD 


In this Section we still consider, as before, a simplified HOD 
which does not distinguish between central and satellite 
galaxies. However, we now consider a cut-off mass for the 
halo richness. Me, which is different from the mass scale Mi 
(In fact, typically Me < Mi — see, e.g.. Tinker et al. (20081). 
Hence, our HOD is: 


Nm = (^ 


MiJ 

Defining the variable: 


(39) 


a o a 

"=r ^2 


and the cut-off: 


M 

Mf 


(40) 


a ( Me 

^ 2 ( AA 


(41) 


the calculations of the previous Section can now be per¬ 
formed in basically the same fashion. The number density 
of halos that can host galaxies is: 


n-h,g — 


M. 


A 

7^ 


/ 0 \ ^0 poo 

i-j j dxx^° [l + (2x)~’’] e~ 


no 


(a) 4r(l + ^o,a:c) + 2-^’r(l-bAo-p,a:c)] , 

(42) 


where, as previously, Ao = —1/2 — 3/(3 -1- n), and we have 
defined no = pmAjMi,i/n. 

A similar calculation leads to the mean number density 
of galaxies: 



X [r(l -|- Al, ic) + 2 ^F(l Xi — p, Xc)] , 


and the galaxy bias becomes: 


bg 


I- + - --TIE-,(441 

= r(l-bAi,a:.)-b2-PF(l-f Al-p,*.)^ ^ 
X [r(2 -b Al, Xc) -f- 2 ^r(2 -|- Al — p, Xc) 

-|-p2“^r(l-b Al - p, a:c)] . 
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Figure 3. Galaxy bias obtained using the Sheth-Tormen for¬ 
malism (blue lines and filled region), and in the Press-Schechter 
formalism (red, long-dashed lines and filled region), as a function 
of the cut-off scale Xc [see Eq. The parameters were allowed 

to range in the intervals 0.9 < q < 1.1, and —2 < n < —1. The 
black (short-dashed) line shows the approximation of Eq. ( |45| l. 


In the limit of Xc ^ 1 the galaxy bias can be considerably 
simplified, in fact: 

lim bg ^ 1 — + 2S~^Xc , (45) 

Xc = OCi 


which is basically the approximate expression we obtained 
in the PS case, where xi played the role of the cut-off mass 
scale. This relationship between the cut-off scale and the bias 
allows us to write, in the limit of high biases, Xc — &c{pg ~ 
l + 5~^)/2. This turns out to be a fairly good approximation 
as can be seen from Fig.[^ Notice that we do not expect the 
model of Eq. (391 to hold for low-biased galaxies, when the 


small-halo mass limit is critical. 

After some algebra, the 1-halo term can be expressed 
in the same way as was done for the PS case: 


T-)! h 

P =q 


rih,, 


(46) 


where: 


9 = 


r(l Ao, Xc) -(- 2 ^r(l -(- Ao — p, Xc) 

[r(i -i- Ai, Xc) + 2“pr(i -|- Ai — p, Xc)] 

X [r(i -(- A2, Xc) 2 ^r(i -(- A2 — p, Xc)^ . 


(47) 


In t he limit of Xc 1 we can use the series expansion of Eq. 
(291 to show that, as in the PS case, the prefactor q —>■ l[^ 
. In this limit reduces to the following expression: 




2 

no 1 - 

a 




(1 + 2 


-Pj.-?) g- 


(48) 


^ In this limit, we see from Eq. ( |46[ | that the 1-halo term of the 
power spectrum inherits a dependence on the number density of 
the halos that contain at least one galaxy. However, the HOD we 
used in this Section and in the previous one make no distinction 
between central and satellite galaxies — in fact, we have simply 
used the typical parametrizations used for satellites. Hence, in 
this context, fih,g should be regarded as the number density of 
halos containing more than one galaxy. 


with Xc depending on the bias in the following way: 

= y(fe9 - 1 + <5c“^)- (49) 

Conversely, in the limit ic <C 1, we obtain the following 
scaling for the one-halo term: 

piH ^ r(i + A2) + 2-^r(i + A2-p) ^ 

^0 (a) ° 2“^r(l -I- Al — p)]^ 

Comparing Eq. ( |36[ ) and Eq. ( |48| l, we see that, for high 
values of the bias, the 1-halo term in the ST model behaves 
in basically the same way as was found for the PS formalism. 


3.3 Semi-analytical model: Tinker mass function 
and realistic HOD 

The analytical approximations of the previous Sections have 
allowed us to obtain simple expressions for the number den¬ 
sities of galaxies and halos, the bias, and the 1-halo terms, 
but we made some strong assumptions — in particular, 
about the simple scaling of halo richness, about the assump¬ 
tion of self-similarity which fixed the scaling of u with mass, 
and about the way in which we cut off the halo richness be¬ 
low some given mass scale. Although the final results may 
have seemed natural and physically sensible, they could have 
been influenced or even driven by these simplifications. 

In this Section we argue that these results are robust. 
We show this by improving the modeling of the previous 
Sections in a number of ways: first, we calculate the mass 
variance a{M) of Eq. 0 from the power spectrum of a 
vanilla-ACDM model; second, we employ t he [Tinker et al)] 
(20081 mass function and halo bias of Eqs. (22l-(23l, which 


are a slightly better fit to the N-body simulations com¬ 
pared to the PS or ST expressions; and third, we use a 
class of HODs which is inspired and calibrated by observa¬ 
tions ( Zheng et al.||20d5 Zheng, Coil &i Zehavi|20OT Zheng 


et al. 20091. We also distinguish between central and satellite 


galaxies — whereas in the preceding Sections we implicitly 
assumed that all galaxies were satellites. 

As for the HOD, we have used the formulas of Eqs. 0- 
([^ for the halo richness of central and satellite galaxies. Typ¬ 
ical values for these parameters are Me — 10^^'® h~^ Mp 
Ml ~ 10^^ ' h-' 
et al.|2009 ). 


Mp, 




In contrast to the previous Sections, where all calcula¬ 
tions could be carried out exactly, here we instead compute 
numerically the galaxy number density of Eq. 0 , the bias 
of Eq. (101, and the 1-halo term of Eq. (12l. In this way 


we can explore basically any point in parameter space, and 
compute the properties of the galaxy models corresponding 
to those points. 

We allow the parameters to vary in the following ranges, 
while keeping always Me < Mi: 


12.75 

< 

login h/Mp) < 14.25 , 

(51) 

13.2 

< 

login -^1 h/Mp < 15.0 , 

(52) 

0.85 

< 

a < 1.15 , 

(53) 

0.85 

< 

<Jg < 1.15 , 

(54) 

1.0 

< 

Kg < 1.3 . 

(55) 


The results for the 1-halo term of the power spectrum are 
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shown in the top panel of Fig. We have split the dif¬ 
ferent HOD models in four groups, according to the num¬ 
ber densities of galaxies. From the bottom, the horizontal 
shaded areas correspond to HODs whose ranges of n~^ fall 


in the intervals 10 


-i+0.15 


> rig Mpc > 10 


for 


L = 3.3, 3.6, 3.9, and 4.2 . The values of for models in 
those four groups are shown as points of the same colors as 
the shaded areas, from lower left to upper right, respectively. 
The dashed black line is the power-law 6g ®. 

We saw in Sections 13.II and 13.21 that the one-halo term 
evolves like a power-law in bg for intermediate values of the 
galaxy bias. This is again what we observe here: the asymp¬ 
totic behavior of the one-halo follows a steep power-law as a 
function of galaxy bias. The growth of the one-halo term is 
constrained by the requirement that Me < Mi, which limits 
the number of galaxies in halos, and imposes an upper limit 
on the bias and the one-halo term. If we relax this physical 
requirement, the power-law evolution of the one-halo term 
continues at higher bias. 

The results obtained in Sections 13.11 and 13.21 were de¬ 
rived under the assumption that all galaxies are of the “satel¬ 
lite” type — an approximation that we did not use here. For 
highly biased tracers the 1-halo term can be comparable to 
Poisson shot noise, which shows that the “central” galaxies 
are less relevant in that limit, providing a motivation for 
that approximation in Sections |3.1| and |3.2| Conversely, the 
fact that drops well below the level of galaxy shot noise 
for low values of the bias shows that, for such types of ob¬ 
jects, the central galaxy plays an important role, as many 
(or most) halos host a single central galaxy. 

In the bottom panel of Fig. we show the results for 
the 1-halo term of the trispectrum of Eq. (14l, in the limit 
A: —>■ 0, fc' —l 0, for the same range of HOD parameters 
used in the top panel. Again, we separate the models in 
groups, according to the number density of galaxies, and 
the horizontal shaded areas denote the different values of 
rig ® for each group. The dashed line is the power-law bg^'^. 

As argued above, the 1-halo term of the trispectrum 
also scales rapidly with bias, approximately as (P^^)®. In 
Fig . [^we show the 1-halo term of the trispectrum [Eq. @1 
against the 1-halo term of the power spectrum. The dashed 
line indicates the scaling (P^^)^'®. 




Figure 4. Top panel: 1-halo term in the HOD of Eqs. §-0- 
The models were split according to the number density of galax¬ 
ies. From left to right, and from the bottom up, the models 
shown have IQ-^ iS > Ug > 10-3'45^ 10-3-45 > fig > IQ-^ '^®, 
10-3.75 > fig > 10-"^ °®, and 10-^ °® > rig > IQ-'^ ®®. In each 
case, the horizontal shaded area marks the corresponding range 
of Ijng. The dashed black line is the power-law feg ®. The stripes 
seen mostly in the lower right corner are just an artifact of the 
grid we used to explore the HOD parameter space. Bottom panel: 
1-halo term of the trispectrum for the same class of HODs. The 
dashed black line is the power-law bg®'®- The horizontal shaded 
areas mark the corresponding ranges of 1 /fig. 


3.4 Simulations 


We further test the results from the previous sections in a 
fully numerical set-up, by using a catalog of halos from the 


DEUS simulations 1^ ( Alimi et al.||2010 Rasera et al.|[2010 


Courtin et al.|20TT I. We use halos detected with the Friends- 
of-Friends algorithm in a box of 648 h~^ Mpc, containing 
1024® dark matter particles, in a ACDM cosmology in agree¬ 


ment with WMAP5 (Komatsu et al. 12009 b We should note 


that the DEUS mass function is very well fit by the Tinker 
mass function, and that the DEUS halo bias is also well fit 
by the Tinker halo bias. 

We populate the halos using the HOD formalism pre¬ 
sented in Eqs. ||^ and ([^, with fixed parameters a. = 0.9, 
cTg = 1.0 and Kg = 1.1. The Me and Mi parameters are 
allowed to vary in the ranges given in Eqs. (|51|l-([5^. 


® http://www.deus-consortium.org 


The number W of central galaxies is either 0 or 1, and 
is drawn from a nearest-integer uniform distribution with 
mean N^. For halos containing a central galaxy, the number 
of satellite galaxies is then drawn from a Poisson distribution 
with mean Ns (since we are taking the A: —f 0 limit for the 
1 -halo term, it is irrelevant where in the halo those satellite 
galaxies are placed). Obviously, halos that do not contain a 
central galaxy have no satellite galaxies. 

With the DEUS halo catalog, and the HOD implemen¬ 
tation described above, expressions such as that for the num¬ 
ber density of galaxies become: 


POO 

J Ml 


dlnM 


duh 

dlnM 


X N{M) 


{iV,[M(/i)] + iV,[M(/i)]} 


^ Nh{M)[Ne{M) + Ns{M)] 


( 56 ) 
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Figure 5. 1-halo term of the trispectrum plotted against the 1- 
halo term of the power spectrum, for the HOD of Eqs. §-0. 
The dashed black line is the power-law The models were 

split according to bias: from lower left to upper right, the different 
colors indicate models with 1.0 < bg < 1.5, 1.5 < 5g < 2.0, 
2.0 <bg < 2.5, and 2.5 < bg < 3.0. 


where V is the volume of the simulation. In order to compute 
galaxy bias, for simplicity we employ the Tinker bias, Eq. 


(231, in Eq. (lOl. 


The resulting scaling between the one-halo term and the 
galaxy bias is shown in the top panel of Fig. The same 
type of scaling found in the previous sections appears here 
— the grey line indicates the power-law cx bg. In order 
to check whether this scaling is sensitive to the precise form 
of the HOD used, we also test a different HOD prescription 
for which: 


M 

Ng oc log -j- for M > Me- 
IVIc 


(57) 


The number of galaxies in each halo is drawn from a Poisson 
distribution with mean given by the formula above. The 
scaling obtained in this case between the one-halo term and 
the galaxy bias is shown on the bottom panel of Fig.[^ The 
grey line now follows 6g'®, a scaling very similar to the one 
found in the previous cases. 

As is the case in the original HOD of the previous sec¬ 
tions, in this alternative HOD prescription the richness is an 
increasing function of halo mass. We also tested HOD mod¬ 
els where the richness becomes constant, or even decreases, 
with increasing halo mass, and in those cases we do not re¬ 
cover the same type of scaling between bias and the one-halo 
term of the power spectrum. However, we do not consider 
this to be a limitation of our study, as we expect realistic 
HODs to display an average number of galaxies monotoni- 
cally increasing with halo mass. 


4 DISCUSSION 

We have shown that, when the galaxy bias is high enough 
{bg > 3), the 1-halo term of the power spectrum grows faster 
as a function of bias than the 2-halo term. We also showed 
that the 1-halo term of the trispectrum grows much faster 
than the 4-halo term. We argue that the 1-halo terms of 
all the N-polyspectra scale faster than the N-halo terms of 
those polyspectra. 



Figure 6. One-halo term and galaxy bias obtained from the 
DEUS simulated halo catalog populated with the HOD of Eqs. 

(top panel), and with the test HOD of Eq. | |57| l (bottom 
panel). In the top panel, the grey line shows the scaling bg, while 
in the bottom panel the scaling is 5g ®. 


We interpret these results in the following way. Galaxy 
bias provides an intuitive physical interpretation of the way 
in which the visible matter distribution is related to the un¬ 
derlying DM distribution. Given a Gaussian field So, taken 
from a distribution whose variance (in Fourier space) is ba¬ 
sically the matter power spectrum, the density contrast of a 
biased tracer on a grid of finite-volume cells can be approx¬ 
imated by a lognormal field, 5g = exp[6g5G — b^aQ/2] — 1, 
where bg is the bias of the tracer, and (Jq is the variance 
of the Gaussian field on the volume of the cell of the grid 
( Coles fc Jones|fl990 l. As the bias increases, the number of 
particles found in the density peaks grow very fast ([Bardeen] 


et al. 19861 — exponentially, in the lognormal model — 


while most of the space is emptied. 

In terms of the Halo Model, we can vary the HOD pa¬ 
rameters in such a way that the mean number density is 
kept hxed while the bias increases. For highly biased galax¬ 
ies we expect to find more objects concentrated in fewer 
halos. This implies that, for very high values of the bias, the 
galaxy 2-point statistics can get a large contribution from 
the statistics of halos hosting two or more galaxies. This 
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argument extends to all the A'^-point statistics: in the high- 
bias limit there will be many galaxies inside the same few 
halos, which means that all the 1-halo terms of the A'^-point 
statistics will become increasingly important. 

In Section 3.1 within the Press & Schechter (1974 ) for¬ 
malism, we used very simple formulas for the halo richness 
and for the scaling of the mass variance (and peak height), 
to show that the 1-halo term of the power spectrum grows 
very fast with bias. A more refined analytical calculation, 
done in Section 3.2 using the Sheth & Tormen (19991 frame¬ 
work, shows basically the same scaling. In Section |3.3| we 
employed the Tinker et al. (20081 mass function, realistic 
HODs, and an exact computation of the mass variance in 
ACDM models, and confirmed that the 1-halo term of the 
power spectrum scales at least as fast as ~ Sim¬ 

ilarly, we showed that the 1-halo term of the trisp ectr um 
scales as fast as ~ b 


A Finally, in Section 


3.4 


used a halo catalog derived from the DEUS simulation, and 
two different kinds of HODs, and again obtained the same 
scaling laws. 

The fact that simple analytical arguments show the cor¬ 
rect scaling of the 1-halo terms with bias is a hint that 
these results should be related to basic properties of Gaus¬ 
sian fields and the matter power spectrum. In fact, trac¬ 
ers of different biases effectively probe different scales in 
the power spectrum: higher/lower masses (and higher/lower 
biases) are related to larger/smaller scales. Since the mat¬ 
ter power spectrum has a power index that ranges between 
n ~ —1 at large scales, to n ~ —2 at small scales, the 
power index which is relevant for the halo masses typical 
of a certain type of galaxy is also an indicator of the bias 
of that galaxy. Conversely, bias also tells us how galaxies 
are distributed among the dark matter halos, and about the 
typical peak height which corresponds to a galaxy with that 
bias. In particular, a higher bias implies that galaxies will 
be more concentrated on fewer halos, enhancing the 1-halo 
terms. Since all these properties can be traced back to the 
slope of the power spectrum, it is not surprising to find that 
different types of galaxies present the same scaling of the 
1 -halo term with bias. 

Since the 1-halo terms are constant on large scales, for 
cosmological surveys that cannot resolve the inner structure 
of halos, these terms enter effectively as additional sources 
of noise and covariance. Conversely, very accurate and com¬ 
plete surveys will be able to detect much better the ampli¬ 
tudes and scale dependences of the 1-halo terms when the 
bias is sufficiently high, implying better constraints on HOD 
parameters. 

In particular, cosmological surveys targeting highly- 
biased objects could be severely impacted by the effects of 
the statistics of counts of the halos hosting those objects. 
Measurements of the power spectrum from these surveys 
should take into account not only the higher effective shot 
noise coming from the 1-halo term of the power spectrum, 
but also the additional contribution to the power spectrum 
covariance coming from the 1-halo term of the trispectrum. 
Similarly, measurements of the bispectrum which employ 
highly biased tracers should ensure that the relevant 1-halo 
terms are properly taken into account. 

Finally, although we worked at 2 = 0, it would be in¬ 
teresting to find out what happens at high redshifts. On the 
one hand, the bias of a given population of tracers is typ¬ 


ically increasing as a function of redshift; but on the other 
hand, linear theory should become a better approximation, 
which means that the 1-halo term should be less important. 
Therefore, at higher redshifts the form of the scaling of the 
1 -halo term should depend sensitively not only on the am¬ 
plitude of the power spectrum, but also on the evolution of 
the HOD parameters. 

The cosmic infrared background (CIB) is a good exam¬ 
ple: in that case (redshifts 2 ~ 1 — 4), the 1-halo term has 
a higher amplitude compared with shot noise, which helped 


constrain the HOD of the sources the CIB (Thacker et al. 


20131. A similar situation may arise as galaxy surveys aim 


at deeper redshifts and higher number densities with ob¬ 
jects such as emission-line galaxies and quasars [e.g., JPAS 
(Benftez et al. 20141; eBOSS (Dawson et al. 20151; DESI 
( Levi et al.|2013 '; PFS ( Ellis et al.|2012 1 ], HI intensity maps 
(with, e.g., SKA^, etc.: the amplitude of the 1-halo terms 
should be kept in check in order to ensure the reliability of 
the forecasts from these surveys. We plan to examine in fu¬ 
ture work how extrapolations of the existing HODs to higher 
redshifts impact the signal-to-noise levels and the cosmolog¬ 
ical constraints for some future galaxy surveys. 
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